Half-metallicity and efficient spin injection in AlN/GaN:Cr (0001) heterostructure 
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First-principles investigations of the structural, electronic and magnetic properties of Cr-doped 
AIN/GaN (0001) heterostructures reveal that Cr segregates into the GaN region, that these interfaces 
retain their important half-metallic character and thus yield efficient (100 %) spin polarized injection 
from a ferromagnetic GaN:Cr electrode through an A1N tunnel barrier - whose height and width 
can be controlled by adjusting the Al concentration in the graded bandgap engineered Ali-^Ga^N 
(0001) layers. 



PACS numbers: 71.20.-b, 73.20.-r, 75.70.Cn, 85.75.-d 

The injection and manipulation of spin polarized carri- 
ers - a key feature determining successful spintronics - is 
highly dependent on phenomena and proccesses occuring 
at the interface between a ferromagnet and a semiconduc- 
tor. It has been shown 1-3 that surface or interface sen- 
sitive approaches must be employed in both theory and 
experiment for the accurate determination of the spin 
polarization in materials proposed for potential magne- 
toelectronic devices. As perhaps the most striking ex- 
ample, Heusler compounds (such as NiMnSb, Co 2 MnGc, 
Co2CrAl) - which possess the appealing half- metallic fer- 
romagnetic behavior in bulk - show significantly reduced 
spin polarization at the surface 2 or, consequently, at the 
interface with a semiconductor 3 . Indeed, first-principles 
calculations have revealed that structural reconstruction 
in the vicinity of the junction essentially alters the elec- 
tronic states and destroys half-metallicity 3 . 

Recently, nitride-based semiconductors, in particular 
GaN and A1N, have attracted increasing attention for 
a number of reasons: (i) The introduction of magnetic 
dopants and achieving ferromagnetism in these materials 
provide complementary functionality to the wide range 
of devices already developed for pure large band-gap 
nitrides 4 , (ii) The shorter bond length and smaller spin- 
orbit coupling in these light-element compounds as com- 
pared to other III-V semiconductors (e.g., GaAs) are pre- 
dicted to give rise to higher Curie temperatures 5 ; indeed, 
recent measurements in Cr-doped GaN and A1N bulk 
showed T c to be over 900 K 6 . (iii) The smaller spin-orbit 
interaction is also one of the main arguments for longer 
(by three orders of magnitude) electron spin lifetimes in 
GaN than in GaAs 7 . (iv) The III-N semiconductors ap- 
pear to be much more resilient to the presence of ex- 
tended structural defects than are other semiconductors 8 . 
While Cr-doped bulk GaN and A1N systems have been 
extensively studied 6,9-11 , both experimental investiga- 
tions and theoretical modeling of magnetically doped III- 
N interfaces - essential for practical spintronic applica- 
tions - are lacking. 

In this Letter, we present results of first-principles 
calculations of Cr-doped AIN/GaN (0001) heterostruc- 
tures focusing on their structural, electronic and mag- 
netic properties. From a comparison of the formation 



energy of the relaxed AIN/GaN heterostructures with dif- 
ferent site locations of Cr, we predict that the magnetic 
impurity segregates into the GaN region which serves as a 
ferromagnet, while A1N is the semiconductor part of the 
interface. Most significantly, we find that the Cr-doped 
nitride-based interface retains the desired half-metallic 
behavior that enables the efficient 100 % spin polarized 
injection across the interface and thus makes this sys- 
tem highly attractive for spintronic applications. Finally, 
we propose a way to control the shape (i.e., height and 
width) of the nitride potential barrier by adjusting the 
concentration x in a Cr-doped interface with a graded 
Alx_ x Ga x N barrier region - a realistic magnetoelectronic 
structure that can be grown epitaxially. 




FIG. 1. Geometry of (2x2) AIN/GaN (0001) heterostruc- 
ture: the large, intermediate and small spheres correspond to 
Ga, Al and N atoms, respectively. Substitutional site loca- 
tions of Cr atoms in the supercell denoted by numbers. 

The AIN/GaN (0001) interface is modeled using (lxl), 
(2x2) and (2x4) supercells of pure or Cr-doped wurtzite 
GaN and A1N (0001) with 3 double layers of each ma- 
terial, cf. Fig. 1. Note that a "double layer" is a 
single nitride layer, i.e., it consists of a group III ele- 
ment and nitrogen, so the (lxl), (2x2) and (2x4) super- 
cells contain 12, 48 and 96 atoms, respectively. For the 
Ali^Ga^N/GaN (0001) heterostructures, we used (2x2) 
supercells with a;=0.00, 0.25, 0.50, 0.75 and 1.00 in the 
layers, i.e., the content of Al (Ga) gradually decreases (in- 
creases) along the (0001) direction until a composition of 
pure GaN is formed. The equilibrium relaxed geometry 
of the structures was determined via total energy and 
atomic forces minimization for the lattice parameters a 
and c and the internal parameter u. Both the GaN and 
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A1N regions of the heterostructure were allowed to re- 
lax with the same in-plane lattice constant 12 , while the 
internal atomic relaxation for each lattice parameter c 
provided different minimized III-N distances along the 
(0001) direction. Note, that during the optimization, all 
atoms were allowed to move in the x, y and z directions. 

TABLE I. Calculated lattice parameters, A, and band 
gap values, eV, for bulk wurtzite A1N and GaN and for the 
AIN/GaN (0001) heterostructure. Experimental data are 
taken from Ref. 16. 





Lattice parameters 




Band gap 






Q, C &exp Cexp 


LDA 


sX-LDA 


Exp. 


A1N 


3.085 4.993 3.112 4.982 


4.34 


6.05 


6.20 


GaN 


3.142 5.194 3.189 5.185 


2.00 


3.35 


3.39 


AIN/GaN 


3.127 15.402 


2.59 


3.98 





We employ the all-electron full-potential linearized 
augmented plane wave (FLAPW) method 13 that has 
no shape approximation for the potential and charge 
density. Cut-offs of the plane- wave basis (16.0 Ry) 
and potential representation (81.0 Ry), and expansion 
in terms of spherical harmonics with I < 8 inside the 
muffin-tin spheres were used. Summations over the Bril- 
louin zone were carried out using 32 special k points in 
the irreducible wedge. In addition to the local density 
approximation (LDA) for the ground state properties, 
we used the self-consistent screened-exchange LDA (sX- 
LDA) method 14,15 , which is known to provide a consid- 
erably improved description of the excited state (optical) 
properties as compared to the LDA or generalized gradi- 
ent approximation (GGA) calculations. 

We first determine the effect of the interface on the 
structural and electronic properties of pure AIN/GaN 
(0001) and compare with those of bulk GaN and A1N. 
Due to the lattice mismatch, both GaN and A1N arc 
found to be mutually strained in the junction, cf. Table 
I: (i) The resulting in-plane lattice constant, 3.127 A, is 
0.5 % smaller (1.3 % larger) than that of the bulk GaN 
(A1N) and is in agreement with the value observed 12 for 
the thinnest AIN/GaN bilayer, 3.134 A. (ii) The renor- 
malized c lattice constant of GaN in the heterostructure 
is larger than that of A1N by 0.195 A which is approx- 
imately equal to the difference between the c parame- 
ters of the bulk materials, cf., Table I. We also found 
that away from the junctions, the Ga-N and Al-N dis- 
tances along (0001) in each sub-unit tend to relax back 
towards their unstrained bulk values and differ by 0.7 
% (0.9 %) from those in the bulk GaN (A1N). Further, 
the structural relaxation at the junctions affects the elec- 
tronic properties, in particular, the optical bandgap of 
AIN/GaN (0001). In Table I, we present the bandgap 
values of bulk GaN and A1N and the (lxl) AIN/GaN 
(0001) calculated using both the LDA and sX-LDA meth- 
ods. As expected, the LDA underestimates the band gap 
of both bulk GaN and A1N by 30-40 %, while the sX-LDA 



gives a considerably improved description - namely, a less 
than 2%-diffcrcncc from the experimental values. For the 
AIN/GaN (0001), we found that the valence band max- 
imum and the conduction band minimum are formed by 
states of the GaN layers; thus, these states determine the 
calculated minimum band gap of the system to be 2.59 
cV (3.98 cV) within the LDA (sX-LDA). In contrast, the 
states of the A1N layer located further away from the 
junctions lie deeper in the valence and conduction bands 
giving a bandgap increase of ~0.4 eV within both the 
LDA and sX-LDA 17 . 

Next, we investigated the effect of the interface on 
the transition-metal impurity location and the magnetic 
properties of the doped structure. To this end, we per- 
formed calculations of Cr doped substitutionally into the 
(2x2) AIN/GaN (0001). (For comparison of the electronic 
and magnetic properties, we also calculated Cr-doped 
bulk GaN and A1N.) The choice of the (2x2) supercell 
is motivated by the recently calculated critical separa- 
tion between two Cr atoms in bulk GaN, ^2.7 A, above 
which ferromagnetic (FM) coupling dominates the anti- 
ferromagnetic (AFM) one 11 . In our supercell case, the 
Cr-Cr distance is 6.3 A. Now, to find the preferred site 
location of the magnetic impurity, we calculated the for- 
mation energies 18 of six relaxed structures with cation- 
substituted Cr, cf. Fig. 1. The results, gathered in Table 
II, allow the following conclusions: 

TABLE II. Comparison of the calculated formation ener- 
gies, eV, and the magnetic moments on the Cr atoms and 
its tetrahedral N neighbors in plane (N p ; anar ) and along 
the (0001) direction (N apica i), fin, for the Cr-doped (2x2) 
AIN/GaN (0001) heterostructures. Site locations of the Cr 
atoms (in parentheses) correspond to those in Fig. 1. 





A1N region 
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Gr(l) 


Cr(2) 


Cr(3) 


Cr(4) 


Cr(5) Cr(6) 


AE, 


+1.825 


+1.817 


+1.827 


0.000 


+0.031 +0.010 


Cr 


2.309 


2.309 


2.308 


2.302 


2.304 2.299 


^planar 


-0.010 


-0.011 


-0.015 


-0.020 


-0.019 -0.016 


N apical 


-0.018 


-0.007 


-0.005 


-0.003 


-0.012 -0.012 



(i) The magnetic impurity segregates into the GaN re- 
gion of the heterostructure which serves as a ferromagnct 
(as shown below), while the A1N region is a non-magnetic 
part of the interface. The difference of the formation en- 
ergies for Cr in the A1N and in the GaN region of the 
heterostructure (~1.8 eV, cf. Table II) shows little de- 
pendence on the distance from the junctions 19 and agrees 
well with that for Cr-doped bulk GaN and A1N (2.0 eV). 

(ii) Although the calculated local magnetic moments 
on the Cr atoms in the AIN/GaN, Table II, are similar 
to those obtained for the corresponding relaxed Cr-doped 
bulk GaN (2.304 fi B ) and A1N (2.313 fi B ), the presence 
of the junctions does affect the magnetic properties of the 
system. This can be seen from a comparison of the in- 
duced magnetic moments on the tetrahedral N neighbors 
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of the Cr atom (i.e., three in-plane N and one apical N 
atoms, cf. Fig. 1), presented in Table II. We found that 
the magnetic moments are larger on the N atoms which 
bond with Ga than on those connected to Al atoms by 
0.006 \ib (or 0.009 fin) for planar (or apical) N atoms, on 
average. This finding can be explained by the fact that 
the states of the A1N layers lie deeper in the valence band 
(as discussed above) and hence they are more screened 
as compared to the GaN layers 20 . Therefore, at the junc- 
tion where the in-plane N atoms bond with Ga and the 
apical N atom bonds with Al, the induced in-plane spin 
polarization is dominant (case 4, Table II). In contrast, 
a higher spin polarization along (0001) - as compared 
to the in-plane magnetic moments - is observed at the 
other junction (case 1, Table II). Accordingly, we found 
the magnetic moments on planar and apical N atoms to 
be similar when Cr is equidistant from the junctions. 




FIG. 2. Geometry of the most energetically favorable Cr 
site location in the AIN/GaN (0001) heterostructure (left) 
and contour plot of the calculated spin density distribution 
within a slice passing through this Cr atom and its planar N 
neighbors (right). Solid (dash) lines in the plot denote pos- 
itive (negative) spin polarization. The calculated magnetic 
moment on the atoms, /^s, is given in parentheses. 

A contour plot of the calculated spin density distribu- 
tion within a slice passing through the Cr atom and its 
planar N neighbors (case 4 in Table II) is presented in 
Fig. 2. In agreement with the itinerant sp-d exchange 
model 21 , the p-orbital of the N atoms which points to- 
ward the Cr atom has its magnetic moment antiparallel 
to the one on the magnetic impurity giving rise to FM 
coupling between two Cr atoms. (Note, that the other 
two N p-orbitals possess positive spin polarization and 
create hopping channels for itinerant electrons.) Indeed, 
from the (2x4) Cr-doped AIN/GaN (0001) supercell cal- 
culations, we confirm strong FM coupling between two Cr 
atoms with an energy difference between FM and AFM 
orderings of 97, 106 and 110 meV/Cr atom for case 4, 5 
and 6, respectively. Accordingly, the effective exchange 
interaction parameters - estimated 22 to be 22, 26 and 28 
meV, respectively - show a similar dependence on the Cr 
site locations, which correlates with the decrease of the 
ratio between the induced magnetic moments on planar 
and apical N atoms, cf. Table II. 

(iii) Regardless of the Cr site location in the GaN 
region of the interface (cf. cases 4-6 in Table II with 



a total energy difference of ~kT), the three Cr-doped 
AIN/GaN systems share similar band structure features 
(the corresponding band offsets are discussed in a sepa- 
rate paper 23 ): The magnetic impurity d states hybridized 
with the p states of the neighboring N atoms form deep 
bands in the nitride bandgap, cf. Fig. 3. For the 
majority-spin channel, partially occupied triply degener- 
ate tig and fully occupied doubly degenerate e g bands are 
located at ~2.5 eV and eV, respectively, above the 
valence band maximum. The exchange interaction splits 
the Cr d states by ~2 eV. These findings are in agreement 
with previous studies 9 ' 11 of Cr-doped bulk GaN. 




FIG. 3. Spin-resolved band structure along the high sym- 
metry directions in the Brillouin zone for the most energet- 
ically favorable Cr-doped AIN/GaN (0001) heterostructure. 
The origin of the energy is taken at the Fermi level. 



Most significantly, we found that the Cr-doped 
AIN/GaN (0001) heterostructures retain the desired half- 
metallic behavior (consistent with the calculated integer 
total magnetic moment of 3 /xs/cell) with a bandgap 
in the spin minority channel of ~2.7 eV, cf. Fig. 3. 
This leads to a complete, i.e., 100 %, spin polarization 
of the conduction electrons and thus makes the system 
attractive for high-efficiency magneto-electronic devices 
- specifically, magnetoresistive tunnel junctions. 

Based on these results, in particular on the preference 
of Cr atoms to substitute Ga rather than Al atoms in the 
AIN/GaN (0001), we model a realistic (in the sense of 
practical applications) magneto-electronic system using 
(2x2) supercells of Cr-doped graded Ali-^Ga^N/GaN 
(0001) heterostructures. From a comparison of the total 
energies for the relaxed structures with different Cr site 
locations, we found that the magnetic impurity prefers 
to substitute the Ga atoms which are located in the pure 
GaN (0001) layers (cf. cases 7-14 in Fig. 4). A signifi- 
cant total energy increase is obtained for the structures 
where Cr substitutes the Ga atoms in the mixed nitride 
(0001) layers (cf. cases 1-6 in Fig. 4). This energy de- 
pendence on the composition suggests the possibility to 
control the width and height of the energy barrier by 
adjusting the Ga/Al ratio in the layers perpendicular to 
the growth direction - an important ingredient for device 
design optimization. 
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FIG. 4. Calculated relative formation energies, K, as a 
function of the cation concentration x in the Cr-doped graded 
Ali-^Ga^N/GaN (0001) heterostructures and geometry of 
the corresponding supercell. The numbers on the graph cor- 
respond to those in the supercell and denote the Ga site lo- 
cations substituted with a Cr atom. 

In summary, based on the first-principles calculations, 
we predict that efficient spin injection can be achieved us- 
ing a ferromagnetic GaN:Cr electrode in conjunction with 
an A1N tunnel barrier. Since the Cr-doped AIN/GaN 
(0001) heterostructures are found to be half-metallic, one 
can expect a pronounced increase of magnetoresistance 
in the nitride-based magnetic tunnel junctions. These 
findings make the system highly attractive for spintronic 
applications and call for experimental investigations. 
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